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Cosmic ray spectral hardening due to dispersion in the source injection spectra 
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Recent cosmic ray (CR) experiments discovered that the CR spectra experience a remarkable 
hardening for rigidity above several hundred GV. We propose that this is caused by the superpo- 
sition of the CR energy spectra of many sources that have a dispersion in the injection spectral 
indices. Adopting similar parameters as those of supernova remnants derived from the Fermi 7-ray 
observations, we can reproduce the observational CR spectra of different species well. This may be 
interpreted as evidence to support the supernova remnant origin of CRs below the knee. We further 
propose that the same mechanism may explain the "ankle" of the ultra high energy CR spectrum. 

PACS numbers: 98.70. Sa,96.50.sb,98.38.Mz 



Introduction — Nearly 100 years after the discovery of 
cosmic rays by V. Hess in 1912, several basic questions on 
CRs, such as the source(s) and the propagation effects, 
are still not well understood. Precise measurements of 
the CR spectra and observations of high energy 7-rays 
and neutrinos, are of great importance to approach the 
answers to these fundamental questions. 

There are several major progresses in the CR measure- 
ments in recent years. The balloon-borne experiment 
Cosmic Ray Energetics And Mass (CREAM) measured 
the energy spectra of the major species from proton to 
iron in the energy range from tens of GeV/nucleon to tens 
of TeV/nucleon with relatively high precision 
remarkable hardening at ~ 200 GeV/nucleon of the spec- 
tra of all species was discovered Q. Most recently the 
satellite experiment Payload for Antimatter Matter Ex- 
ploration and Light-nuclei Astrophysics (PAMELA) re- 
ported the precise measurement about the proton and 
Helium spectra with rigidity from GV to 1.2 TV Q. 
PAMELA data show clearly that the proton and Helium 
spectra deviate from the single power-law function above 
~ 30 GV with a hardening at rigidity ~ 200 GV, which is 
basically consistent with the results of CREAM and the 
previous Advanced Thin Ionization Calorimeter (ATIC- 
2) Q| ■ The hardening of the CR spectra challenges the 
traditional CR acceleration and propagation paradigm. 
Models possibly to explain such a spectral hardening in- 
clude the multi-component sources Q, or the nonlinear 
particle acceleration scenarios where the feedback of CRs 
on the shock is essential (e.g., 0-0] )■ 

Model — In this work we propose that the hardening 
of the observed CR spectra is due to dispersion of the 
injection spectra of the CR sources such as supernova 
remnants (SNRs). A superposition of the spectra of 
many sources with a distribution of the injection spec- 



tra would lead to an asymptotic hardening of the final 
spectra of CRs For example for power-law injec- 
tion spectrum _E~~ 7 , if there is a uniform distribution 
of the 7 index p(j) = _ x 1 t ne total spectrum will 
be J^ 2 £->( 7 )d7 cx (E^ 1 -E~^)/\nE, which will 
asymptotically approach the hardest injection spectrum. 
Such an effect was adopted to explain the "GeV excess" 
of Galactic diffuse 7-rays observed by EGRET [lfl |ll|. 

Observations of X-ray, GeV and TeV 7-rays indicate 
that the SNRs can accelerate particles (electrons and / or 
nuclei) up to very high energies. Those particles are 
thought to be the most probable sources of the Galactic 
CRs. Studies of radio emission spectra of SNRs suggest 
that the spectral indices of accelerated particles have a 
significant dispersion instead of a uniform value 
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The Fermi satellite observed several SNRs in the GeV 
13 - 3] . The 7-ray energy spec- 



band with high precision 
tra and the coincidence with molecular clouds of sev- 
eral SNRs suggest that the 7-rays are most likely of a 
hadronic origin, although the leptonic origin is not ruled 
out. A broken power-law injection of protons with the 
break energy several to tens of GeV seems to describe 
the 7-ray data well [2fJ, |2l| . The fit to 7-rays also indi- 
cates that there is a large dispersion of the accelerated 
CR spectra in SNRs. Assuming the hadronic origin of 
the 7-ray emission of SNRs, the modeling of the GeV- 
TeV 7-ray emission from eight sources detected by Fermi 
and Cherenkov telescopes gives the source particle spec- 
tra 71 as 2.15 ± 0.33 and 72 « 2.54 ± 0.44, for energies 
below and above the break respectively 
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Here we calculate the superposition effect of CR spec- 
tra, assuming that CRs are originated from SNR-like 
sources. The injection spectrum of each source is as- 
sumed to be a broken power-law function of rigidity with 
the break from several to tens of GV. The spectral indices 



2 



are assumed to be Gaussian distributed around some av- 
erage values. The normalization of each source is derived 
assuming a constant total energy of CRs above 1 GeV for 
all sources. Since the particle spectra inferred from the 
7-rays might be different from those leaking into the in- 
terstellar space [22I , the injection parameters are adopted 
through the fit to the CR data instead of the ones inferred 
from 7-ray observations. 



We employ the GALPROP code [23| to calculate the 
propagation of CRs, in the diffusive reacceleration frame. 
The main propagation parameters are Do = 5.8 x 10 28 
cm 2 s _1 . 5 = 0.33, va = 32 km s" 1 and Zh = 4 kpc. 
For each major chemical species we use the superposed 
spectra of many sources as input and calculate its propa- 
gation. The B/C ratio is found well consistent with data 
and insensitive to the source spectrum. 

For rigidity below ~ 30 GV, solar modulation needs to 
be considered. The force-field approximation is adopted 
to model the effect 24J . The modulation potential $ de- 
pends on the solar activity, which varies from ~ 200 MV 
at solar minimum to ~ 1400 MV at solar maximum. For 
the period when PAMELA operates, the modulation po- 
tential was estimated to be 450 — 550 MV According 
to the fit to the low energy observational data (see Fig. 
[T]), we adopt $ = 550 MV for proton and Helium (to fit 
the PAMELA data), and $ = 750 MV for Carbon, Oxy- 
gen and Iron nuclei (to fit the HE AO 3 data). A higher 
modulation potential for HEA03 was also found in [251 ] . 

There is a "knee" of the all-particle spectra at PeV 
energies, which indicates the existence of break or cutoff 
on the CR spectra 0. Such a break or cutoff might 
be due to the acceleration limit of sources, the prop- 
agation/leakage effect from the Galaxy, or interactions 
with background particles [27j . Here we adopt two kinds 
of cutoff/break to model the knee structure of the total 
spectra: a sub-exponential cutoff case with the energy 
spectrum above the injection break i? -72 exp(—R/R c ), 
and a broken power-law case with energy spectrum above 
the injection break i? -72 (1 + i?/ii c ) _1 . In both cases we 
assume that the cutoff/break energy is Z-dependent, i.e., 
the rigidity R c is constant 1 . 

The calculated energy spectra of proton, Helium, Car- 
bon, Oxygen, Iron and the total spectrum, together with 
the observational data are shown in Fig. [1] The param- 



eters of the model are compiled in Table Q] The results 
show good agreement with the data of each major species 
and the all particle one. Here a sub-exponential cutoff 
instead of the standard exponential cutoff is required, in 
order to reconcile the all-particle spectra with the pro- 
ton and Helium data around the knee region. Such a 
sub-exponential cutoff may be originated from stochastic 
acceleration of particles in the turbulent downstream of 
weakly magnetized, collisionlcss shocks [28[. For the cut- 
off case, the expected all-particle spectrum is lower than 
the data at energies above tens of PeV. Such a result is 
consistent with the requirement of a "B component" of 
Galactic CRs [2^]. For the break case, the high energy 
data of the all-particle spectra can be well reproduced 
without the "B component" , as was also shown in the 
"poly-gonato" model [3]. 



TABLE I: Source parameters: injection spectra 71, 72 and 
break rigidity Rb, high energy cutoff rigidity R c and solar 
modulation potential <E>. 

71 72 Rb Rc $ 

(GV) (PV) (GV) 

p 1.95 ±0.20 2.52 ±0.28 [5,30] 0.5 0.55 

cutoff He 1.95 ±0.20 2.50 ±0.33 [5,30] 0.5 0.55 

C,Q,Fe 1.95 ±0.20 2.58 ± 0.35 [5,30] 0.5 0.75 

p 1.95 ±0.20 2.52 ±0.25 [5,30] 0.5 0.55 

break He 1.95 ± 0.20 2.50 ± 0.30 [5, 30] 0.5 0.55 

C,Q,Fe 1.95 ±0.20 2.58 ±0.32 [5,30] 0.5 0.75 

Discussion — In such a simple scenario of dispersion of 
injection spectra of CR sources, the observed harden- 
ing of CR spectra by ATIC, CREAM and PAMELA can 
be reproduced. If the CRs are indeed originated from 
a population of sources instead of a single major source 
(e.g.. 



40|, |4l|), such an asymptotic hardening effect due 



It is reasonable to assume a constant R c for all sources for the 
propagation/leakage models and the interaction models. How- 
ever, the break may suffer from a dispersion in the acceleration 
limit models. We have tested that the result with a dispersion 
of R c can actually be well approximated by the model with a 
proper constant R c . 



to dispersion of source properties is inevitable. The in- 
jection parameters are similar to those inferred from the 
7-ray observations of SNRsi, which might be evidence 
that SNRs are the sources of Galactic CRs below ~PV. 

The injection spectra of different elements are slightly 
different, which might be related to the difference in pro- 
duction of different species [42j. We find that the differ- 
ence between various elements in the injection spectra is 
not as large as that in the observed ones, possibly because 
the interaction strengths of various elements are different 
from each other during propagation. It should be noted 
that the PAMELA data actually showed a sharp break at 
~ 200GV and a gradual softening below the break rigid- 
ity Such detailed features cannot be simply recov- 
ered in the present model, where the gradual hardening 
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FIG. 1: Energy spectra of proton (top-left), Helium (top-middle), Carbon (top-right), Oxygen (bottom-left), Iron (bottom- 
middle) and the all-particle one (bottom-right). The local interstellar spectra are labelled with "LIS" and the observed spectra 
after solar modulation are labelled with "Mod". The solid line in each panel represents a sub-exponential cutoff behavior of 
the high energy spectra around the knee region, while the dashed line is for broken power-law type. References of the data 
are— proton: ATIC-2 0, PAMELA 0, CREAM 0, Tibet 0, KASCADE [13]; Helium: ATIC-2 Q, PAMELA 0, CREAM 
0|, Tibet-BD [H], KASCADE [aH; Carbon: HEA03p3, CREAM 0]; Oxygen: HEA03 [H, CREAM 01; Iron: CRN [3], 
RUNJOB [H, SOKOL [H, JACEE [H, TRACER [lJlIEA03 [3], CREAM 0; all-particle: Tibet-Ill Q. The normalized 
all-particle data are derived by combining all data with a rescale based on the extrapolation of the direct measurements [2r3 [. 



is expected. The same tension also exists for the multi- 
component source model and the non-linear acceleration 
model. We expect that future better measurements of 
wide band spectra by e.g., the Alpha Magnetic Spec- 
trometer (AMS02, [43(1) an d the Larg e High Altitude Air 
Shower Observatory (LHAASO, |44() would help to test 
this model. 

There are some implications of the CR spectral hard- 
ening, for example, the imprint on the secondary parti- 
cles such as positrons [451 ] . diffuse 7-rays, and antipro- 
tons [Zjj ]. Based on our model CR spectra, we calculate 
the predicted hadronic-origin diffuse 7-ray fluxes. The 
total diffuse 7-ray emission consists of hadronic, leptonic 
and the extra-galactic components, which is very compli- 
cated. In the Galactic plane the diffuse 7-ray flux is dom- 
inated by the hadronic component as shown in the con- 
ventional CR propagation models [47[ . The 7-ray yield is 
calculated using the parameterization of pp interactions 
given in fi^ . For the effect of heavy nuclei we employ a 
nuclear enhancement factor cm = 1-84 (4^]. The abso- 
lute fluxes depend on the propagation model and spatial 
distribution of CRs. Here we only discuss the relative 
results. The ratios of the hadronic 7-ray fluxes between 
our model expectation and that of the traditional sin- 
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FIG. 2: Ratio of the hadronic component of the diffuse 7-rays 
between the dispersion scenario expectation and the single 
power-law model. 



gle power-law CR spectrum are shown in Fig. [2] It is 
shown that the 7-ray flux also experiences a hardening 
above ~ 50 GeV. A similar conclusion was also derived 
in [4q |. For positrons and antiprotons we expect similar 
behaviors, although the propagation effect may change 
a little bit the quantitative results. The current Fermi 
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7-ray data and PAMELA antiproton data can not probe 
such a hardening of the secondary particles yet. 

We also note that the dip-cutoff structure is very simi- 
lar to the ankle-GZK structure of ultra high energy CRs 
(UHECRs). Considering the fact that UHECRs should 
also suffer from such a hardening effect if they are from a 
population of sources, we would expect the same mech- 
anism to be responsible for the anklc-GZK structure of 
UHECRs. Following we assume the injection spec- 
trum of UHECRs is a broken power-law function with 
an exponential cutoff. The logarithm of break energy 
\og(Eb/eV) is assumed to be uniformly distributed in 
[17, 1 8] and the cutoff energy is E c ~ 5 x 10 19 eV. The 
spectral index is 2.0±0.2 below and 3.6±0.6 above E\>. 
Assuming a pure proton component we present the 
expected superposed UHECR spectrum in Fig. [3] Wc 
would expect a similar result for any other single chem- 
ical species. Note that the Pierre Auger Observa tory 
data showed a gradual increase in mass composition [53J . 
Since the relative abundance of each chemical species is 
not yet well constrained experimentally, we cannot per- 
form a detailed modeling. In any case, we propose that 
the superposition effect from a population of UHECR 
sources with a dispersion in the injection spectrum pro- 
vides a new ingredient to model the UHECR spectrum. 
Since the interactions between UHECRs and the back- 
ground photons are unavoidable if UHECRs are produced 
at cosmological distances [5ol | , we may in turn expect that 
UHECRs are produced locally or even in the Galaxy 
if the mechanism proposed in this work is responsible for 
the shape of the UHECR spectrum. More generally, it is 
possible that both the superposition and the interaction 
effects are in operation to give the ankle-GZK structure 
of UHECRs. 
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FIG. 3: Calculated energy spectra of UHECRs for pure pro- 
tons, compared with the HiRes data [53 |. 
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